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The microscopic mechanism of the metal-insulator transition is studied by orbital-resolved ^*V NMR spec¬ 
troscopy in a prototype of the quasi-one-dimensional system VjOo. We uncover that the transition involves a 
site-selective d orbital order lifting twofold orbital degeneracy in one of the two VOe chains. The other chain 
leaves paramagnetic moments on the singly occupied d^y orbital across the transition. The two chains respec¬ 
tively stabilize an orbital-assisted spin-Peierls state and an antiferromagnetic long-range order in the ground 
state. The site-selective Mott transition may be a source of the anomalous metal and the Mott-Peierls duality. 

PACS numbers: 71.30.+h, 76.60.-k, 75.50.-y, 75.40.Gb 


The metal-insulator transition (MIT) involves several in¬ 
stabilities in quasi-one-dimensional (quasi-ID) systems 111]. 
A well-known example is a tetragonal (rutile R) phase VO 2 
in which a competition of Mott and Peierls instabilities has 
been actively debated for decades ffl. A low-temperature 
nonmagnetic insulating (monoclinic M\) phase has a strong 
dimerization with d orbital ordering along the chain, sup¬ 
porting the orbital-assisted Peierls mechanism of the MIT 
| 10(] . The complex phase diagram including a partially para¬ 
magnetic insulating (monoclinic M?) phase under the slight 
V oxidation lEl or uniaxial stress ifl^ shows degenerated 
metastable states arising from the competition of electron cor¬ 
relation and electron-lattice coupling. The underlying physics 
remains a fundamental issue for strongly correlated systems 
and practical applications to the Mott transistor El. 

A Wadsley phase VeOn is another quasi-lD system QB- 
in which the microscopic mechanism of the MIT has 
been debated for half a century. It has been also extensively 
studied as a potential cathode in a Li-ion battery IBB]. 
There are two frustrated zigzag chains including three vana¬ 
dium sites, V(l)-V(l) and V(2)-V(3) chains along the b axis, 
as shown in Fig. [1] A first-order MIT occurs at Tmi = 150 K, 
accompanied by a drop of the magnetic susceptibility 
[T]l and a weak structural dimerization along the chain ylH. 
The lattice symmetry is lowered from a high-temperature 
monoclinic C2/m phase to a low-temperature monoclinic Pc 
phase, doubling vanadium sites [H . Whereas the nominal va¬ 
lence is the V-O bond valence sum suggests 

nonuniform valence states V(i)+4-26, V(2)+'* ®, and V(3)'*’"* ®' 
even in the metallic phase H 18|. Similar to the M 2 phase of 
VO 2 EI remains paramagnetic in the insulating phase, im- 
plyin g a sp iri-singlet formation for only a part of the vanadium 
sites 12311241] . Residual local moments exhibit antiferromag¬ 
netic long-range ordering at Tw = 50 K Eilal- 

It is noted that the metallic phase of VgOis exhibits sev¬ 
eral anomalous features, x strongly depends on temperature 
T above I’m: (Fig. [T]|, implying the presence of local mo¬ 
ments. The photoemission spectroscopy shows the depressed 
density of states at the Fermi leyel ll26l429ll . suggesting that 
the spectral weights are partially shifted into lower Hubbard 
bands due to the strongly renormalized electron correlation. 
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FIG. 1. Magnetic susceptibility x measured with a SQUID magne¬ 
tometer for a single crystal of VeOo {Ho = 1 T || and b'). The lower 
and upper insets show crystal structures in the metallic (monoclinic 
C2/ra) and insulating (monoclinic Pc) phases, respectively 0]. 


These features are in sharp contrast to those of VO 2 showing 
the Pauli paramagnetic susceptibility and the appreciable den¬ 
sity of states at the Fermi level El , although Tmi in VeOia is 
suppressed to the lower temperature compared to Tmi = 342 
K in VO 2 , possibly due to competing ground states. 

In this Letter, we address the microscopic mechanism of 
MIT and the origin of the anomalous metal in the quasi-ID 
vanadate. We focus on the orbital degrees of freedom as a key 
to understanding the MIT for the d electron system, since the 
electron correlation can be renormalized into d orbital occu¬ 
pations IB] . However, the experimental probe of orbitals has 
been limited in simple transition-metal compounds Eol] . and 
the direct observation of the d orbital for several atomic sites 
remains challenging. Here, we explore NMR spectroscopy 
as a site-selective orbital probe by utilizing anisotropic hy- 
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perfine interactions. We show a dramatic d orbital ordering 
across Tmi for one of the vanadium sites, while remaining as 
paramagnetic d^y character for the other sites. 

NMR measurements were conducted on a sin gle c rys- 
tal of VeOn made by chemical transport using TeCU lllhll and 
coated with epoxy to prevent crushing at the structural transi¬ 
tion. Frequency-swept NMR spectra were obtained from spin- 
echo signals using a njl - t - njl pulse sequence {njl - 1 
//s, T = 5 - 30 ps) with a 0.3 MHz step at a constant mag¬ 
netic field Hq - 9.401 T. For a given orientation of the crystal 
in Hq, each vanadium site V(0 (/ = 1, 2, 3 in the metallic 
phase and i - la, lb, 2a, 2b, 3a, 3b in the insulating phase) 
gives a set of the NMR spectra at resonance frequencies 
= [1 (n = -3, -2,..., 2,3), where 

are defined as relative shifts from vq = 105.23 MHz and 
is the nuclear quadrupole splitting frequency for the nuclear 
spin 7 = 7/2. The maximum \6v\ gives the nuclear quadrupole 
frequency vq - 3e^qQfh2I(2I - 1) with Planck’s constant h, 
the nuclear quadrupole moment eQ, and the electric field gra¬ 
dient (EFG) eq - Vzz- The nuclear spin-lattice relaxation rate 
Tj' was obtained from a single exponential fit to the recovery 
curve after an inversion n pulse at a magic angle {6v - 0) for 
each V site. 

Figure 13 a) shows the NMR spectra for Hq || c* (a'*). 

In the metallic phase above 150 K, three sets of the spectrum 
with the characteristic K and 6v are reasonably assigned to 
three vanadium sites in reference to the expected valences and 
the symmetry axes of VOe distortions. Spike lines located 
around vq (blue) come from V(2) close to W^^(3d^) 1^ . The 
others with the negative K and small (red) or large (green) 6v 
[Figs. ISa), (b)] are respectively from V"^+-like V(l) or V(3) 
with a weak or strong VOe distortion. Below Tmi, the number 
of spectra increases due to the site doubling. A difference in 
the V(3a) and V(3b) spectra indicates a slight charge transfer 
between the vanadium sites across Tmi. The most noticeable 
change occurs for the V(l) spectrum; Two sets of the sharp 
quadrupole split (colored with gray and red) appear with the 
sizable 6v and vanishing K, indicating the d orbital order, as 
shown below. 

First we focus on the local susceptibility for each vana¬ 
dium site based on K measurements (Fig. SI, Supplemental 
Material 1^ 1. obtained from K is shown in Fig. [3 
exhibits a strong thermal variation in the metallic phase. It 
shows a discontinuous increase (decrease) for V(3a) [V(3b)] 
at Tmi. The T dependence behaves similar to that of the spin 
S - 1/2 ID Heisenberg antiferromagnet (TY = 2 7S, ■ Sj) 
with an exchange coupling 7 = 77 K for;t''^*^^^ and be¬ 
low 140 K and with 7 = 69 K for;if''^^^ above 150 K. In con¬ 
trast, exhibits rather weak T dependence in the metallic 
phase and becomes suppressed in the insulating phase, leav¬ 
ing the small Van Vleck susceptibility [< 5 x 10“"^ emu/V(l) 
mol]. The result gives clear evidence for site-selective spin- 
singlet formations. Namely, the V(l) chain exhibits the MIT 
from a (Pauli) paramagnetic metal to a nonmagnetic insulator, 
as seen in VO 2 d ■ becomes negligible below Tmi: 

The on-site d spin is absent on V(2a/2b), rendering V(la/lb) 




FIG. 2. (a) V NMR spectra of VgOo [77o = 9.401 T || c*(fl'*) for 
the metallic (insulating) phase], measured with t = 6 and 10 - 30 
ps for the lower and higher frequency parts, respectively. The red, 
blue, green, and gray lines are assigned to V(l) [V(la)], V(2) [V(2a), 
V(2b)], V(3) [V(3a), V(3b)], and V(lb) based on the angle depen¬ 
dent profiles (Fig. 4). (b) K and (c) 6v obtained from the V NMR 
spectra. 


and V(3a/3b) effectively half-filling. 

The dynamical property measured with 7/^ is another man¬ 
ifestation of the site-dependent correlation [Fig. [3b)]- In 
contrast to weakly cotTelated metals with the T-linear 1/Ti, 
IjTiT for V(3) increases upon cooling and then becomes fur¬ 
ther enhanced below Tmi for V(3a) and V(3b) due to the evo¬ 
lution of antiferromagnetic fluctuations toward T^. For V(l), 
IjTiT exhibits a large drop at Tmi, consistent with the be¬ 
havior of;^''^b) Reflecting the dilute spin density, IfTiT for 
V(2) is an order of magnitude smaller than those of the other 
sites. Thus, the site dependence in the spin correlation be¬ 
comes highlighted below Tmi. 

To give deeper insights into the microscopic origin of the 
MIT, the predominant d orbital occupation is investigated by 
measuring the angular dependence of K and Sv (Fig. |4]l. 
Since the magnetic susceptibility is isotropic in vanadates, 
the K and 6v anisotropies are governed by the anisotropic 
hyperfine couplings expressed with the equivalent operators 
q ,7 = \{LiLj + LjLi) - 6ijL^ (i,j - x,y,z) lUlit], where L rep¬ 
resents the total orbital angular momentum in the local coordi¬ 
nates (x,y,z) corresponding to (a, b, c*) or (c', b\ a'*) (Supple- 
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FIG. 3. (a) Site-dependent local susceptibility obtained from the 
isotropic Knight shifts and (b) IjTiT in VeOis. Solid curves 
are fitting results of the ID Heisenberg model, whereas the dashed 
one is an average of and^''^^'’*. 


mental Material 03611 1. One can elucidate orbital occupations 
by extracting q,j. For a singly occupied dxy, dy^, or d^x orbital, 
the components of q,j exhibit the cyclic angular dependence 
[Fig. Ufa)] with a minimum along the symmetry axis of the 
orbital. For metals, q,j represents the electron cloud distri¬ 
bution for the Hund-coupled bands crossing the Fermi level. 
With increasing electron correlation, the spin-polarized lower 
Hubbard band may signihcantly contribute to the spin suscep¬ 
tibility i^ . On the other hand, q,j obtained from 6v measures 
the net d occupation below the Fermi level. 


As for the metallic phase, has a minimum for the c* 
axis [Fig. Efb)], in good agreement with the dxy anisotropy, 
whereas exhibits out-of-phase angular dependence with 
a maximum for the c* axis and the depressed anisotropy. q,y is 
extracted from the spin parts of the axial anisotropy 
and the in-plane asymmetry expressed with a linear 

combination of three t 2 g orbital components dxy : t/yj : d^x - 
a \ p \ j {a + p + y = 1), where Kl-^ - (2K* - - Ky)l2>, 

Kin - andK®^ = H-Kp/3 (Supplemen¬ 

tal Material Obll l. The obtained occupation ratio is listed in 
Table H] It turns out that V(l) has an admixture of the dy^^ and 
d^x components. It is consistent with the depressed value of 
vq -0.1 MHz due to the degenerated orbital occupation (Ta¬ 


ble S3, Supplemental Material 0611 1. The results are also in 
qualitative agreement with a density functional theory (DFT) 
calculation (a : p : j - 0.0/0.23/0.38) |34|. In clear contrast, 
the V(3) site is mostly polarized to dxy, suggesting strong elec¬ 
tron correlations under the crystal field. 


TABLE I. Axial anisotropy and in-plane asymmetry 

^In/^iso’ nuclear quadrupole frequency vq, the EFG asymme¬ 
try T], and the orbital occupation in VjOo, where the occupations are 
obtained from the K anisotropy for V(i) (i = 1, 3, la, lb, 3a, and 
3b) and from the nuclear quadrupole frequency for V(la) and V(lb) 
(Supplemental Material (3^). 


site 


kL/kl 

vq(MHz) 

7 

dxy . dyr . d-x 

V(l) 

-0.38 

0.21 

0.11 

0.40 

0.10:0.56:0.28 

V(3) 

0.77 

0.03 

0.46 

0.26 

0.77:0.13:0.10 

V(la) 

- 

- 

0.49 

0.47 

0.11:0.18:0.71 

V(lb) 

- 

- 

0.64 

0.20 

0.13:0.78:0.09 

V(3a) 

0.70 

-0.05 

0.44 

0.69 

0.74:0.16:0.10 

V(3b) 

0.74 

-0.04 

0.54 

0.76 

0.76:0.14:0.10 


On the insulating phase, the angular profiles of K and 6v 
for V(2a/2b) and V(3a/3b) are similar to those of the metallic 
phase. It confirms that the orbital states of these sites have 
less impact on the MIT. The most remarkable change across 
I'm! is the emergence of the huge vq sites, 0.49 MHz for V(la) 
and 0.68 MHz for V(lb) 1^ . pointing to d orbital ordering 
on the V(la)-V(lb) chain. The minimum of 6v appears along 
the b' (c') axis for V(la) [V(lb)], in agreement with the pre¬ 
dominant d^x (dyz) occupation [Figs. 4 (a), 4(e)]. Namely, the 
dyzidzx degeneracy in the metallic phase is completely lifted 
by the orbital order. The Sv anisotropy gives the largely polar¬ 
ized occupation for V(la) and V(lb) (Tables U and and Sup¬ 
plemental Material Table S3 after subtrac ting the lattice 
contribution (Table S2, Supplementa Material 13611 1. We thus 
conclude that the MIT in the present system is attributed to 
site-selective orbital ordering accompanied by the spin-singlet 
formation. 

The obtained predominant orbitals are depicted in Fig. |5] 
Owing to the dilute spin density on V(2), the dxy orbitals on 
the V(3) sites are weakly coupled each other on the V(2)-V(3) 
zigzag chain and carry the local moments even in the metal¬ 
lic phase. The degenerated dy^jd^x orbitals on V(l) form the 
conducting chain, which are lifted into the singly occupied dy^ 
or dzx state below Tmi. A weak dimerization along the V(la)- 
V(lb) chain produces the spin-singlet pairs. The n bonding 
between orthogonal d orbitals should be rather weak. This 
feature is in clear contrast to the d orbital order in VO 2 where 
the direct d-d cr-bond stabilizes strong spin-singlet pairing 
M- To explain the vanishing spin susceptibility with a siz¬ 
able spin gap (A > 150 K) in VeOn, one may have to con¬ 
sider significant superexchange interactions through the O Pz 
orbital. Such interactions help to lift the^/^Z^ degeneracy, 
analogous to the Jahn-Teller mechanism 13411 . In context, the 
MIT involves the spin-Peierls transition assisted by orbital or¬ 
dering and electron-phonon coupling. The above observation 
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FIG. 4. (a) Calculated q,; proportional to the on-site dipole hyperfine coupling and EFG for dy., and orbitals. The orthogonal coordinate 
axes (x,y,z) are taken on the VOe octahedron and parallel to (a, b, c*) or (c', b', a'") in Fig. [T] (b) - (e) Angular dependence of the V Knight 
shift K and the nuclear quadrupole splitting frequency 6v at 170 and 100 K in VsOb. The curves are fitted results with the sinusoidal functions 
of Volkov’s formula 
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FIG. 5. Configurations of the predominant orbital occupations in 
the (a) metallic and (b) insulating phases of VeO^. The bold lines 
indicate spin-singlet pairs. The arrows on V(3a)/V(3b) denote aniti- 
ferromagnetically ordered spins below 50 K 1^ . 


of the paramagnetic local moment on V(3a/3b) provides evi¬ 
dence for the Mott insulating ground state of V 6 O 13 with an¬ 
tiferromagnetic ordering i24ll. whereas a theory beyond the 
DFT +U band calculation 1341] is required to describe param¬ 
agnetic Mott insulators. It is remarkable that the orbital oc¬ 
cupation and the exchange interaction do not show significant 
changes for V(3) across keeping the Mott insulating fea¬ 
tures. The results showing the partially localized vanadium 
sites are consistent with the depressed spectral weight at the 
Fermi level and the significant weights well below the Fermi 
level as observed in the XPS measurements 11271] . 


earth nickelates based on a DFT - 1 - dynamical mean field the¬ 
ory calculation 1371] . while it has not been confirmed exper¬ 
imentally. Our observations of the site-dependent d occu¬ 
pations and spin susceptibilities give microscopic evidence 
for the site-selective Mott localization. Many transition- 
metal compounds with inequivalent atomic sites have site- 
dependent local environments and a potential to the site- 
selective localization, giving anomalous bad metals with the 
enhanced resistivity and magnetic susceptibility. Even under 
the weak crystal field, the emergence of the large orbital polar¬ 
ization in the metallic phase may feature the evolution of the 
lower Hubbard bands. The orbital-resolved NMR provides 
the site-selective local probe of the orbital order due to strong 
mass renormalization and can be applied to extensive materi¬ 
als 


To conclude, we have presented the observation of the site- 
selective metal-insulator transition, which accounts for the 
anomalous features in the historical quasi-one-dimensional 
correlated system VgOia. The orbital polarization and local 
spin susceptibility are elucidated for three vanadium sites with 
different electron correlations. We uncover that the metal- 
insulator transition is dominated by orbital ordering and spin- 
singlet pairing for the itinerant vanadium site with the orbital 
degeneracy, whereas the magnetic site with the singly occu¬ 
pied orbital remains with paramagnetic moments. Thus the 
spin-Peierls state and antiferromagnetic order coexist in the 
insulating ground state. The results demonstrate the compet¬ 
ing instabilities and the reasonable ways of symmetry break¬ 
ing in the frustrated one-dimensional system. 


The site-selective Mott transition has been proposed in rare- 
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Supplemental Material for Site-Selective Mott 
Transition in a Quasi-One-Dimensional 
Vanadate VgOn 

In Supplemental Material, we describe the procedure to ob¬ 
tain the site-dependent local susceptibility and the d orbital 
occupation from the NMR measurements in V 6 O 13 . 


A. LOCAL SUSCEPTIBILITY 


The net magnetic susceptibility (emu/V-mol) observed by 
the bulk measurement is expressed as an average of the local 
susceptibilities [emu/V(0-mol], namely, 

A (SI) 

i 

for the metallic phase {i - 1, 2, 3), and 

(S2) 

i 


for the insulating phase (/ = la, lb, 2a, 2b, 3a, 3b) of V 6 O 13 . 

consists of the spin and orbital components, 

Aorb*’ for d electrons and the diamagnetic one for core elec¬ 
trons ^dia, 


V(i) V(0 , V(0 , 

A =As +Aorb +Tdia- 


(S3) 


For 3d compounds, Xorh is governed by the temperature- 
independent Van-Vleck paramagnetic susceptibility ;t'vv- We 
redefine after subtracting ;^fdia - -6 x 10“® emu/V-mol. 

NMR probes as the isotropic Knight shift Ki^o consist¬ 
ing of the spin part 



(S4) 

and the orbital part 


= A”Vorb/NA/tB 

(S5) 


with the spin and orbital hyperfine coupling constants. A® and 
A”*, the Avogadro number A/ 4 , and the Bohr magneton jUb- 
For vanadates. A® is usually negative due to the predominant 
core polarization interaction. 

We measured the Knight shift for the three orthogonal di¬ 
rections, (a, b, c*) in the metallic phase, corresponding to (c', 
b', a'*) in the insulating phase, as displayed in Fig. SI. Since 
X would be isotropic for vanadates due to weak spin-orbit cou¬ 
pling, the anisotropy of K arises from the anisotropic dipole 
hyperfine coupling. The isotropic part, the average of the three 
orthogonal components, Ki^o = (Ka + K/, + Kc.)/3, gives 
by determining A®, as mentioned below. 

In the insulating phase of go to 

zero at low temperatures. Then;^^ consists and 

^V(3a/3b)^ Where ;r''<3a/3b) ^ ^v(3a) ^^V(3b) ^^V(3a) +;y(3b)^/2. 

From the Knight shifts at 60 K, = 0.03% and = 

® ISO ISO 

0.31%, we obtained = 4 x 10“® emu/V(la)-mol and 
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FIG. SI. Temperature dependence of K obtained from the V NMR 
spectra for Hq parallel to the three crystal axes in VgOo. 


;^.yOb) ^ 4 x 10-5 emu/Vnb)-mol by using A°* = 2 (r'^) = 
461 kOe/jUB for V^^ (Ref. qI). The residual 4 /''( 3 a/ 3 b) plotted 
against K^'^^VSb) ^ [^v( 3 a) + ^v(3b)]/2 in Fig. [SSa). The 
linearity yields the hyperfine coupling constant A® = -9.5 
T/jUb. The crossing point with Eq. (S5) (the line with the 
positive slope) gives K°* - 0.4% and = 5 x 10“^ 

emu/V(3a/3b)-mol. It is natural to use the common A® and 
in V(3a) and V(3b) with the similar orbital states, giv¬ 
ing 4 /V( 3 a) fromK''(3^> and [Fig. 3(a)]. By 

subtracting and from 4 ^ below 140 K, we confirm 
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thatscales to 
with a reasonable A® = -12 T//iB, as shown in Fig. S2(b). 
Using the common A® for V(la) and V(lb), ^v(ib) 

respectively giveand;if'*'^'*’^ [Fig. 3(a)]. 

As for the metallic phase, three vanadium sites have nega¬ 
tive Knight shifts, reflecting the predominant Using the 
common A® and Xoib to the insulating phase, we obtain 
and from and Finally, the subtraction of 

and;^^''® from;i^ gives 

B. ANISOTROPIC KNIGHT SHIFTS 

The d orbital occupation is reflected as the anisotropic 
hyperflne coupling. The magnetic part of the nuclear spin 
Hamiltonian “TTmag is written as 

■^^mag = -P[kS ■ I -h ■ q ■ I] - jNfitl ■ a - I, (S6) 

where S and I are electron and nuclear spin operators, re¬ 
spectively, ft - hjln is the Planck’s constant, ^ and 

P represents the product of the Bohr magneton the nu¬ 
clear gyromagnetic ratio yN, and the ionic radial average fac¬ 
tor (r“^) = 3.684 a.u. for (3£/') ||T1]. The first term in 


Eq. (S 6 ) represents the core polarization interaction giving 
the negative isotropic Knight shift proportional to the spin sus¬ 
ceptibility with a coefficient k. The second term arises from 
the dipole hyperflne interaction proportional to both;^fs and the 
orbital polarization q. The components q,y {i,i - x,y,z) are 
expressed with the equivalent operators 

q,y = ^(LiLj + LjU) - 5ij\? (S7) 

with the total orbital angular momentum L. The third term 
corresponds to the temperature-independent Van-Vleck term 
arising from a cross term of L and the Zeeman interac¬ 
tion under magnetic held with the components - 
4ju| (r^^) Afiv, where A^v is the second-order mixing proba¬ 
bility of the ground and excited states. 

In a 3d^ case, the principal components of q, y are given by 
q^z — '^^xx — “7qyy — 6 for dxy, ^xx — T ^Qzz — ^ ^ttr 

dy^, and q^^, = -2q^-;c = - 6 for d^^ IliH]. For arbitrary 

occupations, q, j can be written as a linear combination of three 
orbital components, (q^^^, q^.^,, q^j.) = 3(-a + 2(5 - y,-a - p + 
2y, 2a-jS-y), where : dy^ : d„ - a : p : yia+p + y - 1). 
In VgOn, the local coordinates (x, y, z) are taken parallel to the 
crystal axes (a,b,c*) or {c',b',a'*) defined in Fig. 1, , which 
are close to the V-O bond directions. For metals, q,/ reflects 
the d orbital anisotropy for the electrons at the Fermi level. 

By measuring the angular dependence of K (Fig. 4) we 
obtained the principal components (Kx, Ky, Kz), as displayed 
in Table S1. We And that the principal axes are directed close 
to the crystal axes for V(0 (/ = 1, 3, la, lb, 3a, and 3b), in 
good agreement with the V-O bond directions {x,y,z). We 
can evaluate from the K-x plots and extract for each 
axis. Then the Knight shift is decomposed into the isotropic 
part K? ^ = {K^ -H K® -H K^)I3, the axial part K*,; = ~^x~ 

Ky)l3, and the in-plane asymmetric part Kl^ - (K^ - Ky)/2, 
which are expressed as = -kPxs, Klx - j^Xsi'^a-p-y), 
and = ^Pxsi-P + y), respectively, for x = X, y = Y, and 
z - Z. One can cancel out Xs and extract the occupations 
by dividing Kl^ and by Using the typical k - 0.5 
for vanadates II, we obtain = 7 ( 2 a - p - y) and 

KlJKl^ = ^i-P+y). TheresultsforV(l), V(3),andV(3a/3b) 
are summarized in Table I. 

C. NUCLEAR QUADRUPOLE FREQUENCY 

The nuclear quadrupole splitting Sv in presence of the elec¬ 
tric held gradient (EFG) manifests the spatial d electron distri¬ 
bution via the electric hyperflne interaction. The electric hy¬ 
perflne interaction "Kei between the nucleus (the nuclear 
spin 7 = 7/2, the quadrupole moment Q - -0.05 x 10“^"^ cm^) 
and 3d electrons is expressed as|^ 

r\ 2 /~\ 

Hi - ^( 2 ^_l) g(U^)[I-q-I]. (S 8 ) 

Here q has the same form as in the magnetic dipole hyperflne 
interaction in Eq. (S 6 ) and reflects the anisotropic distribution 
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where qzz represents the principal components of the orbital 
polarization given by 3(-a + 2j3-y) for V(la), 3(-a-j3 + 2y) 
for V(lb), and 3(2a-/3-y) for V(l), V(3), V(3a), and V(3b). 

The experimental results of 6v measured for the three crys¬ 
tal axes are shown in Fig. The sign of dv is selected so 
as to satisfy 6va + dvb -i- - 0 and agrees with the angu¬ 

lar dependence of K. The maximum |dv| gives the nuclear 
quadrupole frequency vq proportional to the sum of the EFG 
from the outer ions and the d electrons. 


ve = (i-roo)ve‘ + (i--R)ve. (SiO) 

The outer ion contribution is enhanced by the core elec¬ 
trons with the Sternheimer antishielding factor, jaa, whereas 
the d electron contribution is reduced with a shielding fac¬ 
tor, R ~ 0.35 for vanadates Q]. 

The outer contribution to EFG, V™' (m, n - x,y, z), is ex¬ 
pressed as 


yOUt 

mn 




,3m;n; - 


-)■ 


(Sll) 


with the charge on the /-th ion, qi, and the distance between 
the V sites and the i-th ion, r,. We evaluated V™' with a point- 
charge approximation within a 200A sphere using the crystal 
structure at 290 K for the metallic phase 13] and at 90 K for 
the insulating phase |3] . Here we postulated the valence state 
of V4+ for V(l) [V(la), V(lb)] and V(3) [V(3a), V(3b)], and 
V5+ for V(2) [V(2a), V(2b)]. Then the lattice contribution to 
the quadrupole splitting frequency, dv”'"', can be calculated as 


3eyout 
c out _ * mn 

“ h2i{2i-\y 


(SI 2) 


The results are shown in Table S2. We can evaluate joo as ~ 
-5.5 by comparing the averaged vq ~ 240 kHz (Table S3) and 
V™' ~ 37 kHz (Table S2) for V(2a) and V(2b) without on-site 
d electrons. Using the common jco for other vanadium sites, 
we extract for V(l) and V(3) with the on-site d electron, as 
listed in Table S3. From Eq. (S9), the orbital occupations are 
elucidated (Table S3). The results are in qualitative agreement 
with those obtained from the Knight shift in Table I. 


FIG. S3. Temperature dependence of Sv measured for Hq parallel to 
the crystal axes in VeOo. 


of 3d electron clouds. The nuclear quadrupole frequency Vg 
due to 3d electrons is given by 



e^Q 

lhl{21 - 1) 


{r ^)C|zz. 


(S9) 
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r(K) 

site 

Kx{%) 

Ky(%) 

Kz(%) 

exO 

0x(°) 

ez(°) 


170 

V(l) 

-1.9 

-2.3 

-0.9 

-5 

11 

-84 

92 


V(2) 

0.05 

0.05 

-0.60 

47 

98 

-13 

1 


V(3) 

-1.1 

-1.2 

-3.4 

-2 

-92 

-72 

6 

100 

V(la) 

0.06 

-0.28 

0.06 

75 

10 

87 

80 


V(lb) 

0.55 

0.15 

0.21 

-74 

23 

-10 

80 


V(3a) 

-1.0 

-0.9 

-2.6 

0 

80 

40 

0 


V(3b) 

-1.5 

-1.4 

-3.8 

0 

80 

40 

0 


TABLE I. Principal components of the Knight shift, {Kx,Ky,Kz), and the directions of the X and Z axes with respect to the crystal axes 
(a, b, c") or (c', b',a"‘) at 170 and 100 K in VgOn, where 6 is measured from c* and 0 from a in the ab plane. 


r(K) site (kHz) Svy (kHz) ^Vz (kHz) gx(°) Ml fe(°) Ml 


290 

V(l) 

V(2) 

V(3) 

10.2 

-18.4 

-23.5 

12.1 

-19.3 

-33.0 

-22.3 

37.7 

56.6 

90 

90 

90 

0 

90 

90 

0 

4 

8 

0 

1 

0 

90 

V(la) 

-11.1 

-24.8 

36.0 

78 

-12 

-12 

-3 


V(lb) 

-13.9 

-20.1 

34.1 

92 

79 

79 

-12 


V(2a) 

-18.1 

-20.7 

38.9 

75 

-2 

-15 

0 


V(2b) 

-13.8 

-22.4 

36.3 

73 

-4 

-17 

2 


V(3a) 

-26.5 

-31.1 

57.7 

91 

87 

-20 

2 


V(3b) 

-24.8 

-31.2 

56.1 

90 

88 

-20 

1 


TABLE 11. Nuclear quadrupole splitting frequency calculated in a point charge approximation at 290 and 90 K in VsOo, and the directions of 
the X and Z axes. The X direction is taken perpendicular to the Z axis giving vq. 


r(K) 

site 

5vx (kHz) 6v 

r (kHz) dy 

z (kHz) 6x0 M°) SzC) Ml 


rf a 


7 

170 

V(l) 

33 

78 

-111 

-5 

11 

-84 

92 

183 

1.0 0.24 0.33 

0.43 


V(2) 

163 

84 

-247 

78 

80 

-2 

43 

- 

- 

- 

- 


V(3) 

404 

56 

-460 

0 

85 

18 

5 

827 

0.32 0.77 

0.18 

0.04 

100 

V(la) 

130 

362 

-491 

75 

10 

87 

80 

725 

0.18 0.11 

0.18 

0.71 


V(lb) 

252 

382 

-636 

-74 

23 

-10 

-80 

857 

0.10 0.13 

0.78 

0.09 


V(2a) 

96 

-187 

283 

58 

92 

-28 

29 

- 

- 

- 

- 


V(2b) 

69 

-129 

198 

58 

92 

-28 

29 

- 

- 

- 

- 


V(3a) 

60 

340 

-440 

0 

80 

40 

10 

815 

0.76 0.77 

0.19 0.04 


V(3b) 

68 

480 

-540 

0 

80 

40 

10 

904 0.81 0.81 

0.19 

0.0 


TABLE III. Principal components of the nuclear quadrupole splitting frequency, the directions of the principal axes, v^, t;'*, and the orbital 
components for the metallic (170 K) and insulating (100 K) phases in V^Ob, where Vq is defined as the maximum absolute value and the 
in-plane EFG anisotropy rj = \6vx - Svy\Ivq. X and Y are taken as the directions perpendicular to the Z axis giving vq. 



